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1. Introduction 



In recent years an exciting area of theoretical physics has been unveiled through the discovery 
of the AdS/ CFT correspondence and its numerous generalizations, [1]. Holographic techniques 
have proven a powerful analytic tool in studying the qualitative properties of strongly interacting 
physical systems with applications ranging from describing the physics of strongly interacting quark 
gluon plasmas to modeling collective condensed matter phenomena such as superconductivity and 
superfluidity. 

Despite the success of the AdS/CFT correspondence in studying holographic gauge theories, namely 
field theories with known dual supergravity backgrounds, there are many realistic theories of great 
phenomenological importance, such as QCD, which do not have an explicit holographic dual. This 
is the main reason of the limited direct quantitative applications of the correspondence. However, it 
turns out that under extreme external control parameters, such as: temperature, chemical potential 
and external electromagnetic fields, different gauge theories exhibit similar properties. Therefore it 
is natural to apply holographic techniques to study phenomena which are known to be of universal 
nature. 

One such phenomenon is the magnetic catalysis of mass generation. At weak coupling it has been 
extensively studied using the conventional perturbative field theoretic techniques [2] . A holographic 
study of this effect has been performed in [3] , where the case of flavored A/" = 4 Yang-Mills theory has 
been explored. Further holographic studies of this phenomenon have been addressed by numerous 
authors ^ (for a comprehensive review of the literature look at ref. [4]). At present all such studies 
are in the limit of the so called "quenched" approximation when the number of fundamental fields 
Nf is much smaller than the number of color degrees of freedom N^,. On the gravity side this 
corresponds to the probe limit for the flavor D7-branes. 

In this paper we undertake first steps towards a holographic description of the phenomenon of 
magnetic catalysis of mass generation beyond the "quenched" approximation. This suggests the 
construction of a supergravity background with a fully backreacted set of flavor branes coupled to 
a non- vanishing S-field. This construction amounts to a non-supersymmetric background. 
Ideally such a background would correspond to a set of fully localized flavor branes. The gravity dual 
of a fleld theory with unquenched flavor is coming through the solution of the equations of motion 
with brane sources. It is the presence of these sources which typically modify the Bianchi identities 
and through their contribution to the energy-momentum tensor, the Einstein equations also. If 
the flavor branes are localized, the sources contain Dirac delta functions and, as a consequence, 
solving the equations of motion is, in general, a difficult task. In the context of the AdS/CFT 
correspondence, the search for localized solutions was initiated in [6], where the D3-D7 intersections 
were discussed, while a lot of progress was reported in subsequent years, (see e.g. [7]). 
We could overcome the technical difficulties of localized flavor brane embeddings by considering 
conflgurations which are partially smeared ^ along appropriately chosen compact directions. Super- 
symmetric backgrounds corresponding to partially smeared Karch-Katz like embeddings, [9], have 
been constructed in [10]. In general these backgrounds possess a hollow cavity in the bulk of the 

^For a concise review in the case of holographic duals to the Dp/Dq-brane intersections look at ref. [5]. 
^For a detailed review on the smearing approach see [8]. 
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geometry where the supergravity solution is sourced solely by the color branes. The radius of this 
cavity is related to the bare mass of the fundamental flavors [11]. In the limit of vanishing bare 
mass the cavity shrinks to the radius of the compact part of the geometry and the supergravity 
background has an essential singularity at the origin of the non-compact part of the geometry. In 
both cases the dilaton field diverges at large radial distances. This corresponds to the Landau pole 
that the dual field theory develops in the UV, due to its positive beta function /3(A) oc Nf/Nc- 
One can imagine that a non-supersymmetric background interpolating between the supersymmetric 
backgrounds corresponding to massless flavors in the UV and massive flavors in IR would describe 
a dynamical mass generation. The radius of the hollow cavity would correspond to the dynamically 
generated constituent mass of the fundamental flavors. 

A promising framework for the construction of such a geometry has been developed in [12]. In 
this paper, a ten dimensional black-hole solution dual to the non conformal plasma oi J\f — A 
Yang-Mills coupled to Nf ^ 1 massless flavors has been presented ^ . The D7 flavor branes are 
smeared homogeneously and extend along the radial direction up to the black hole horizon. The 
smearing procedure reduces the flavor symmetry group from U{Nf) to f/(l)^-^ and allows for a 
simple way to account for the backreaction of the flavor branes. This in turn allows one to explore 
the "unquenched" regime of the dual fleld theory. In the zero temperature limit, the resulting 
backgrounds coincide with those found in [15].^ More precisely, the authors of [12] consider the 
smearing of a general non-supesymmetric "fiducial" embedding of the flavor brane for the purpose 
of obtaining a perturbative non-extremal black hole solution. However it turns out that obtaining 
even a perturbative solution in the general case is technically too difficult and the authors obtain 
a perturbative finite temperature solution for the particular case of massless flavors (when the 
"flducial" embedding is trivial) . 

In this paper we will construct a perturbative non-supersymmetric backgroTind with a non- vanishing 
B- field, which corresponds to an external magnetic field coupled to the fundamental degrees of 
freedom of the dual gauge theory. Following the approach of [12] we will consider the case when 
the "fiducial" embedding is trivial. Note that the results of the probe hmit case considered in [3], 
show that this embedding is unstable since it corresponds to vanishing constituent mass of the 
fundamental flavors. In the probe limit this instability is present for any non-vanishing value of 
the external magnetic field. Therefore one may expect that the background obtained by smearing 
such unstable embeddings would also be unstable. However the effect of backreaction is that the 
theory develops a Landau pole. This suggests the existence of an extra energy scale in addition 
to the energy scale associated to the external magnetic field. Therefore we may expect that the 
supergravity background would be characterized by a non-zero critical value of the 5-ficld (Her) 
below which the background is stable. Our studies on the stability of the background using a probe 
D7-brane confirm this expectation. 
Let us summarize the content of our work: 

In Section 2 we introduce useful notation and the supergravity ansatz that is appropriate for the 
description of the backreacted supergravity background discussed above. We proceed by construct- 

■'All the hydrodynamic transport coefficients of tfie modei were anafyzed in [13], while the addition of a finite 
baryon density was presented in [14] 

"^Other solutions employing the smearing technique appear in [16]. 
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ing an effective one dimensional action and obtaining tlic supergravity equations of motion. After 
that we solve those equations perturbatively to first order an appropriate small parameter, namely 
e* oc A=f We are able to present our solution in a closed form. In the limit of vanishing magnetic 
field our solution reduces to the supersymmetric background constructed in [15]. 

In section 3 of this paper we introduce a probe flavor brane to the supersymmetric background 
obtained in[15]. We first consider the case of supersymmetric embeddings and perform a holographic 
regularization of the probe brane action. Our findings enable us to propose an AdS/CFT dictionary 
apphcable at the finite UV cut off of the theory. Next we consider D7-brane embeddings with fixed 
worldvolume gauge field. On the field theory side this corresponds to a constant external magnetic 
field coupled only to the fundamental degrees of freedom introduced by the probe brane. Using the 
proposed AdS/CFT dictionary we obtain numerical plot of the fundamental condensate versus the 
bare mass parameter and explore the effect that the presence of backreacted massless fiavors have 
on the effect of magnetic catalysis of mass generation. Our findings suggest that the effect is to 
make the magnetic catalysis less efficient. 

In Section 4 we generalize the study of Section 3 for an external magnetic field coupling to both the 
probe and the backreacted fundamental fiavors. To this end, we consider a probe D7-brane in the 
non-supersymmetric background constructed in Section 2. We identify the value of the non-trivial 
S-field at the finite UV cut off as the value of the external magnetic field in the corresponding 
field theory. Our study of the classical embeddings of the probe brane suggest that for sufficiently 
low values of the perturbative parameter e^, < Scr the stable phase of the theory exhibit dynamical 
mass generation. However for sufficiently large values of = ecr the theory is unstable. This 
instability is manifest as a diverging slope of the plot of the parameter of the IR separation (related 
to the constituent mass of the probe fundamental fiavors) versus the parameter e*. We find further 
evidence for this instability by studying the meson spectrum of the theory and identifying appro- 
priate tachyonic modes. We interpret the instability of the probe as refiecting an instability of the 
supergravity background. 

In the Conclusion section of the paper we discuss briefiy our present results and their possible 
extensions. 



2. Constructing the background 

The main goal we want to address with the gravity background of the following subsection is 
the study of the phenomenon of mass generation in magnetic catalysis. The field theories we are 

interested in are realized on the intersection between a set of Nc color D3-branes and a set of Nf, 
homogeneously smeared, flavor D7-branes. This is mainly the construction that was studied in [12] 
and in order to accommodate the phenomenon of mass generation we will impose an additional 
coupling between the fundamental fields and an external magnetic field. 

The color D3-branes are placed at the tip of a Calabi-Yau (CY) cone over a Sasaki-Einstein manifold 
X^, where the latter can be expressed as a C/(l) fiber bundle over a four dimensional Kahler-Einstein 
base (KE). In the absence of fiavor branes the color ones source a background whose near horizon 
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limit is the AdS^ x X5 and the dual gauge theories arc superconformal quivers. ^ 
The flavor D7-branes introduce fundamental matter in the dual field theory. They extend along 
the radial direction of the background, wrap a submanifold X3 of X5 and at the same time smear 
homogeneously over the transverse space [18, 19]. The smeared distribution is taken in such a way 
that the isometrics of the fibcred Kahler-Einstcin space are kept unbroken and allows to write an 
ansatz where all the unknown functions just depend on a single radial coordinate. The D7-brane 
embedding is described by a constant profile, implementing massless flavor flelds in the dual gauge 
theory. As a general feature of all the D3-D7 setups, the dilaton runs and blows up at a certain 
radial distance, corresponding to a UV Landau pole in the dual gauge theory [15]. 
The presence of both kind of branes will deform the ten dimensional space-time, a product of a 
four dimensional Minkowski space with a six dimensional CY cone, through a (self dual) F5 and an 
Fi RR flelds. The coupling of the flavor D7-branes with the external magnetic held will be realized 
by the simultaneous presence of a B2 NS held and its electric dual C2 RR fleld, along the gauge 
theory directions. 

In the next subsection we will be precise about the mathematical expression of every one of the 
forms that we described above. 



2.1 Ansatz 

Our notation will follow closely [12, 14, 15]. The action for the of Type IIB supergravity coupled 
to Nf D7-branes, in the Einstein frame, is given by the following expression 

S^SjjB + Sfi, (2.1) 

where the terms of the Sub action are 



«~ 23! 2 5! 



1 1 ^-$172 

2 3!' ^(-^J 



9k2 



C4 A i/3 A F3 , (2.2) 



and the action for the D7-branes takes the usual DBI+WZ form 

Nf 



J d^xe^^- det(G' + e-*/2jr) - j Cq A {e'^)s_^ 



(2.3) 



with jF = B + 27ra'F. In those expressions B denotes a non-constant magnetic fleld, F the world- 
volume gauge field and the hat refers to the pullback of the quantities, along the worldvolume 
directions of the D7-brane. The gravitational constant and D7-brane tension, in terms of string 
parameters, are 

1 TV 

9s 



9k2 



(27r)7^2Q,'4 



(2.4) 



■'For X5 = the CY is the six dimensional Euchdean space and the dual field theory is Af = 4 SYM 
For = T^'^ the CY is the singular conifold and the dual theory is the Klebanov-Witten quiver [17] 
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The ansatz for the metric that we wiU adopt is inspired by [12] and has the following form 

dslo = h-'^ [-dt^ + dxl + b{dxl + dxD] + [b^S^F^da^ + S^dslp2 + F^{dr + Acp^)"^] , (2.5) 
where the CP' metric is given by 

dsQp2 — \dx^ + J cos' ^{d9^ + sin' 9d(f^) + ^ cos' ^ sin' ^{dip + cos 9d(p)'^ & 

Acp2 = - cos^ —{dip + cos 9 d(p) . (2.6) 

The range of the angles is < (%, 6*) < tt, < 93, r < 27r, < ip < iir. The ansatz for the NS and 
the RR field strengths will be the following 

B2^H dx^ Adx^ , JdtA dx^ , 

F^^Q.il + *)e{S^), Fi^Qf{dT + Acp2), F3 = + A Fi (2.7) 

where ^(S's) is the volume element of the internal space and QcQf are proportional to the number 
of colors and flavors 

(2^Ys.a- ^ VoHX,)g. ' 

In our case X5 = and the volume of the three sphere is 27r'. The fact that the flavors are 
massless is encoded in the independence of on a, (see [10, 15]). All the functions that appear 
in the ansatz, h,b, S, F,^, J H, depend on the radial variable a. In the convention we follow, 
S SzF have dimensions of length, b,h,JhH are dimensionless and a has dimension length"^. The 
ansatz for the F3 RR field strength is determined by the ansatz for Fi (see Appendix A for more 
details). Finally the function b in the ansatz for the metric reflects the breaking the of 5*0(1,3) 
Lorentz symmetry down to 5*0(1, 1) x 50(2). 

The equations of motion and Bianchi identities following from (2.1) will be provided in full detail 
in Appendix A. 

2.2 Effective action and the equations of motion 

Since all the functions that participate in the ansatz for the solution, (2.5) & (2.7), depend only on 
(7, it is feasible to describe the system in terms of a one-dimensional effective action. In order to 
achieve that we insert all the ingredients of the ansatz in (2.1) and integrate out all the variables 
except (7 arriving to 

^e// = ^yA.cia (2.9) 
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where Vi,3 is the (infinite) volume of the Minkowski space and Cid is given by the following expression 



1 /h'V fS'V F'S' 



Producing (2.10) we have not made use of the WZ term, since it does not depend on the metric 
or the dilaton. Its effect has been taken into account through the expression for (see [15, 12]). 
The precise expression of the WZ term is needed in producing the equations of motion explicitly 
from (2.1) and it is this point that the smearing procedure enters the field and imposes the following 
replacement 

... I nA... , (2.11) 

Ma J Mio 



Nf 

where is a form orthogonal to the D7-branes and we call it smearing form. The mathematical 
formula for the smearing form is given through the Bianchi identity for F(i) and it is 

dFr = -gs^2. (2.12) 

We will provide more details on the smearing of the DBI part of the action in Appendix A, where 
we will present all the details about the equations of motion and Bianchi identities of (2.1). Since 
the potential J enters the effective action only through its derivative, it corresponds to a "constant 
of motion". This new parameter is related to the value of the magnetic field close to the boundary 
through the equations of motion for F3, coming from the lOd supergravity. We will fix this constant 
of motion in the following way 

^ ^ -QM. ^ J' ^ {H - H^) . (2.13) 

The next step is to use equation (2.13) to eliminate J', in favor of H^, from (2.10) after performing 
the following Legendre transformation 



(2.14) 



The Euler-Lagrange equations will be calculated from the new, transformed, action (2.14). Defining 
the following auxiliary (dimensionless) expressions 
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we can write the equations of motion in the following compact way 



(2.16) 



62 2Qf hS^F^ 1 



8^Mo,h) ^ -Qi- - 



-(2.17) 



a^(log S) = -2 b^F^S^ + 6 b^F^S^ - 



QfcH^F^S'^ 1 



+ =-e^H^Q}hS'(3s 



(2.18) 



a^(log F) = 4 fe^F^^^ - ^ (1 + Q/ e'*6'5' + 



2„8 QfH^hS^F^ 1 e-*/ii/'2/32 



+ 



62 



(2.19) 



Q)e^^b^S^ + 



AQfb'^e'^S^F^ 



1 

2 P 



(2.20) 



e-^hir_ 



e^QlHhS^ + QcJ' + 



, AQfHhS^F^ 



(2.21) 



It is straightforward to check that the above set of equations, together with (2.13), solve the full set 
of Einstein equations, provided the following "zero-energy" constraint is also satisfied 



b' fh' S' F' 



+ 



(2.22) 



+ AQfb''e^F^S''pi-^^''-^ 



le-'^hH"^2 1 



62 ■ 2 

This constraint can be thought of as the aa component of the Einstein equations. Differentiating 
(2.22) and using (2.13) & (2.16)-(2.21) we get zero, meaning that the system is not overdetermined. 

2.3 Perturbative solution 

The system (2.13) & (2.16)-(2.21) allows for a perturbative solution along the lines of [14]. There, 
the dimensionless parameter e*, which is connected to the position that the dilaton blows up, was 
introduced and subsequently used as an expansion parameter. In order for the solution to be valid 
on a large energy scale this parameter has to be a small number and in principle the smaller the 
value of the number the larger the range of the energy scale. After defining A* as the 't Hooft 
coupling at the energy scale r* the parameter e* is expressed in terms of the physical quantities as 



Qfe' 



Vol{X^) ^ Nf 
A* 



(2.23) 
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and particularly for our case that X5 = we have = g^A*^. The parameter can be thought 
of as a flavor-loop counting parameter in the dual field theory. 

Obtaining a perturbative solution, in terms of of the system (2.13) & (2.16)-(2.21), we need 
to impose the following two requirements to fix the constants of integration in the order by order 
expansion of the solution. The first one is that the geometries should coincide with the extremal 
ones in the absence of the magnetic field and the second is that the functions F, 5 & $ should 
correspond to the expressions given in [14], at the energy scale r = r^,. We will redefine the radial 
variable a in such a way that the warp factor keeps the standard AdS form 



h 



r-4 



& 



(2.24) 



therefore to first order in we have the following expression 



4^4 



72 H 



ar 



at 



+ 



96 H 



a. 



a. 



192 



{al - 1)^ log 



ar + 1 
a^ — i 



- - 1) log 



ar, + 1 



ar 



17 



144^4 



{ar - arj 



16 H 



{ar - 1) log 



+ 1 



ar 



(«r, - 1) log 



ar. + 1 



ar 



+ 



ar 



144 H 



(2.25) 



where another set of auxiliary dimensionless functions are in order 



ar=\l + 



4^4 



& = W 1 + 



Art 



(2.26) 



Having (2.25) to transform the solution of the system (2.13) & (2.16)-(2.21) from the one holographic 
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coordinate to the other we arrive to the following expressions for $, H, J' and b 



{ar + l)(ar. - 1) 
{ar - l){ar, + 1) 



(2.27) 



H — 



1 - 



a, 



a^hl 
a2-l 



— a. 



- 1 



Ctr — 1 



(2.28) 



+ 



az 



log 







^4 - 




* . 





J' 



2r 



ar 



Oil + 1 



OL. 



2 +1 ^4 



1 



— (Xr 



0^1 - 1 



+ 



a. 



2 

2-1 



log 



ar + 1 
aj. — 1 



log 





+ 1" 


,^8 




- 1_ 


^8 



(2.29) 



6=1 + 



a,. — 1 



Integrating (2.22) we have for F & 5 

e 



r 
16 



a, 



("r - 1) - "r. (a', - 1] 



9rf 



r 


^4 


r 


+ 8 




+ 4 . 



Or 



r ar 



+ 



32 



(a^ - l)^log 



a, 
+ 1 

CKr — 1 



2-1 rf a? _ 1 



ar 



while decoupling (2.18) & (2.19) 



F-S ^ 



12 



a, 



+ 



5 r 



+ 



^2 "^r* (o^r* 



log 



1) 



a, 



3/2 



a: 



r 

' ■4- 



(2.30) 



(2.31) 



(2.32) 
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2.4 Validity of the perturbative solution 

The perturbative solution we present in the preceding section needs to be supplemented with a 
hierarchy of scales. In terms of the radial coordinate r (or a) our solution has two pathological 
regions that we will describe in this section both qualitatively & quantitatively. 
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As explained in the introduction, the background we have constructed consists of a sea of massless 
flavors. This means that the flavor branes are stretched down to the bottom of the geometry and 
the charge density is highly peaked at r = 0. This is the reason behind the curvature singularity 
that appears in the origin.^ Another way to understand this pathology is through the restoration 
of the U{Nf) symmetry in the deep IR and the subsequent increase of the effective string coupling 
in such a way that the smearing approach breaks down. Careful choice of the integration constants 
when obtaining the perturbative solution together with the presence of the magnetic field, produce 
a less severe IR singularity but by no means cure it. 

A crucial point is that due to the infrared divergency of the background the perturbative solution 



presented in equations (2.25^ 
particular the Jacobean 



da 
dr 



(2.32) is not valid below certain radial distance above the origin. In 

is vanishing at some r//j(e^,). As one may expect if we increase the 

number of backreacted flavors the value for the radius ^/^(e) increases. One can imagine that higher 
order corrections in could reduce the value r/i?(e), however the intrinsic infrared divergency of 
the background at the origin suggests that riR remains non- vanishing in any perturbative solution. 
In Sections 3 and 4 wc will use this radius as a parameter characterizing the infrared applicability 
of our probe-branc analysis. 

There are two ways of either avoiding or hiding the infrared singularity. The first one is to pull 
the flavor branes away from the origin, while keeping the radial symmetry, by introducing massive 
flavors. Then the distance between the color and the flavor branes is interpreted as a mass for the 
fundamentals. The second possibility is the addition of temperature to the background, which has 
effect of hiding the singularity behind the horizon, see [12]. 

A common characteristic of all the imquenched holographic backgrounds constructed so far is the 
appearance of a point in the deep UV that the dilaton diverges. This behavior is a signal for a 
Landau pole in the dual gauge theory.^ The perturbative solution we constructed in the previous 
section is valid until r — r^, where r* denotes an (arbitrary) UV cutoff scale ^ ~ ^uv, well below 
the energy scale of the Landau pole. Additionally the energy scale of should be well above the 
IR scale, in order for the UV completion to have only negligible effects in the IR physics. The 
main characteristic of this scale is that the UV details of the theory do not affect the IR physical 
predictions. This feature is reflected by the independence of all physical quantities (up to suppressed 
contributions) on the position of r* but only on IR parameters. The quantitative outcome of the 
above analysis is that e^, ^ 1, [12, 14], and as we increase the number of flavors the effective physical 
region for the solution between the IR & UV energy scales decreases. 



3. Probing the supersymmetric background 

In this section we introduce a probe D7-brane to the supersymmetric background of [15], which 
is the limit of vanishing magnetic fleld for the perturbative supergravity background constructed 
above. We consider both supersymmetric and non supersymmetric embeddings. The latter have a 

very nice visualization of this construction can be found in [8]. 
^The first counter example is the addition of flavor D6-branes to the AB JM [20] , where the smeared unquenched 
supergravity solution has a good UV behavior. 
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fixed value of tlie U (l)-gauge field strength along the coordinates X2 and 2:3, namely F23 = const. On 
the field theory side this corresponds to an external magnetic field coupling only to the fundamental 
degrees of freedom introduced by the probe D7-brane. We will contrast those results to the ones of 
Section 4, where the external magnetic field couples to all the fundamental degrees of freedom. Let 
us begin by studying the properties of the supersymmetric embeddings. 

3.1 Holographic renormalization 

In this subsection we propose a holographic renormalization scheme similar to the one reported in 
[21] for asymptotically AdS gravitational backgrounds. Our goal is to come up with a prescription 
to calculate vacuum expectation values (like the fundamental condensate) for the supergravity 
background of [15], which suffers from UV divergencies. In the limit of vanishing magnetic field the 
first order perturbative solution is given by [12] 



+ e, log - + C (e^) , 



24 



1 + 



3rf 



6 = 0, h = 
+ 0{el) & So 



(3.1) 



1 + e, 



24 



Ir^ 
3^ 



+ O {el) . 



The lagrangian of the D7-brane probe is given by 



5* 



where J\f = 8T7 Vol(S'^) Vi,3 with Vi,3 the volume spanned by (t, a;*) and Vol(S'^) = 27r^ the volume 
of the unit wrapped by the probe brane. In the next step we expand the classical embedding 
around according to % = Xo(^) + e*Xi(^) aiid solve perturbatively the EOM obtained from (3.2). 
The zeroth order lagrangian is given by 



= Vr3cos^^fl + 4'o(r)^l 
Af 8 2 V 4 V 



1/2 



(3.3) 



which is just the lagrangian for a D7-brane probe in pure AdSs x space time analyzed in [9] 
(in slightly different coordinates). The solution for a supersymmetric embedding corresponding to 
flavor with bare mass m/{2'Ka') is given by 



m 



XoW) = 2 arcsin 
r 

The solution for Xi(r) corresponding to the first order correction to (3.4) is 

n [r^-m^ + 12r4log(^)] 



(3.4) 



36 r^Vr 



(3.5) 
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In the next step we substitute (3.4) and (3.5) in the lagrangian (3.2), obtain the "on-shell" lagrangian 
Cci to first order in e* 



^cl 1 / 2 2\ I ^* ^ 



2 2 4 1 fTn\ 

12 m r log y — j 



(3.6) 



+ 4 r^- 15 rV^ + m^ (r^ + 3r^) -36rMlog(^^^ 



and then integrate it from m to r*. Note that since Xi("^) = the classical embedding of the 
D7-brane, to first order, closes at r = m above the origin. It is convenient to define 



S,i = J drCa = Si^^ + e^Sil^ + O {el) , 

m 

and after that expand in powers r* 



(3.7) 



o(0) 


1 . 2 

32^^* 


— m^) 




4 

^* 


2 2 




~ 288 ^ 


576 



32 



2 2 

rtm 



m 

16~ ^ 32 ' 



5 + 12 log 



m 
192 



1 + 4 log 



O 



(3.8) 
(3.9) 



The divergent and finite terms can be canceled [21] by the addition of the following counter terms 
at r = r* 



(3.10) 



where 7^1, and are appropriately chosen coefficients and 7 is the determinant of the induced 
metric on the r = const slice. Note that we have integrated along all the compact directions (their 
contribution is included in M) and we have added a factor of e* to all the counter terms. To first 
order in e* we have the following divergent and finite contributions from the counter terms 





^4 
' * 








Arim 


#26*- 




Lf 


16m'^ 







m 



3i?4 *9i^4 

4 



rl + 



m 



1 + 12 log 



5 m 
9^ 



1 + 12 log T- 



(3.11) 
(3.12) 
(3.13) 
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One can check that the foUowing set of coefficients 



#2 



34 
5R^ 



1 



(3.14) 
(3.15) 
(3.16) 



cancels all the divergent and finite terms in the first order of the "on-shell" action 5"^°^ + e^S^l\ 
An interesting observation coming from (3.14), (3.15) and (3.16) is that the coefficient of all the 
counter terms has the same dependence, namely (1 + e*/9). In fact one can check that this 
remains true also for the second order corrections, suggesting that resummation is possible. Indeed, 
it turns out that if one considers the non-perturbative form of the solution the dependence of the 
"on-shell" action factorizes. Let us briefiy sketch the derivation. The supersymmetric background 
is determined by the following set of first order differential equations [12] 



S = a'-^e" 



1 + e* 



6 



P* - P) 



F^a'-^e" [1 + e,(p, - p)]-^ 
$ = -log[l + e*(p* -p)] , 



1 + e* ( - + p* - p 



dh 
dp 



1 + ( - + - p 



(3.17) 

(3.18) 
(3.19) 

(3.20) 



where a new radial coordinate p has been introduced. The relation between r and the new radial 
variable p can be obtained as an expansion in powers scries of e* and to first order is given by [12] 



72 



[e^p-^p* _ 1 + i2(p, - p)] 



+ O (el) 



(3.21) 



Note that for = one has simply r = a'^e^, while this relation also holds at r = r* for arbitrary e* 
Therefore at = the supersymmetric embedding x(p) corresponding to Xo('") from (3.4) is 



x(p) = 2 arcsin 



(3.22) 



The parameter eP'' is related to the bare mass of the fundamental degrees of freedom introduced by 
the probe. It turns out that the supersymmetric embedding at finite e* is still given by equation 



14 



(3.22)^, see [12]. The next step is to evaluate the "on-shell" action Sd and for this we note the zeroth 
order term in (3.8) which suggests the following -S^^^ /{Na'"^) = e^*{e'^P* - 6^^5)^/32. The classical 
action to all orders in e* has the following expression^ 




(3.23) 



Remarkably the dependence of the "on-shell" action factorizes when it is written in the radial 
coordinate p. Furthermore, since at r = the relation between the two radial coordinates is 
= a' 2 6''* the dependence of the "on-shell" action remains unchanged if one defines mo = 
a'^ef" — m -\- 0(e*) as a bare mass parameter. This suggests that to all orders in e* the counter 
terms are given by equation (3.10) with the following coefficients #i, #2, #3 

#. = -^^(l + -..j , (3.25) 

The above considerations suggest that one can consistently regularize the "on-shell" action of the 
probe D7-brane provided one keeps the parameter e* sufficiently small to ensure that the finite 
cut off r* = a'^eP* is sufficiently far from the Landau pole of the theory and at the same time is 
sufficiently large so that terms of order 0{e~''*) can be ignored in calculating the regularized "on- 
shell" action. Therefore we propose that the usual AdS/CFT dictionary can be applied at r = r*, 
namely one can expand 

sm -^—^ = \- — + ( — 

2 n n V* 

and identify uiq = a'^e^i as the bare mass parameter and c oc (qq) as the fundamental condensate. 
Note that for the supersymmetric embedding of (3.22) we have c = and thus the fundamental 
condensate vanish (as it should for supersymmetric theory). In the next subsection we will consider 
non-supersymmetric embeddings with fixed U{1) gauge field strength F23 = const and apply the 
proposed dictionary to study the fundamental condensate of the theory. 

3.2 D7-brane probe with fixed ?7(l)-gauge field. 

In this subsection we introduce a D7-brane probe with fixed U{1) gauge field strength F23 = 
Ho/2Tin' . On the field theory side this corresponds to coupling the J\f = 2 hypermultiplet with 

^Note that equation (3.21) suggests that m = r{pq) = a'^e^" + 0(e). 

^We refer the reader to the appendix of the paper for more details on this calculation. 



(3.27) 
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a constant external magnetic field HQ/2'Ka' along the Xi direction. The effective action of the probe 
brane in the Einstein frame is 



5 = -T7 y" d^xe^^-det(g + e-^/^B2) +Tt j P[C(8)] , (3.28) 

where C(8) is the background Ramond-Ramond form sourced by the smeared flavor branes. One 
can show that the effective lagrangian is 



= -e COS -(^cos - + -sm -j (^1 + J (l + e H,h) 

' ^U'^S^'cos"^ , (3.29) 



32 2 



where the functions F hh are given in equations (3.17)-(3.20). Note that in (3.20) the equation 
of motion for h is in an integral form but it can be obtained in terms of incomplete gamma functions. 
We refer the reader to the appendix for the technical details. 

It turns out that for the numerical analysis it is convenient to define the following dimensionless 
quantities 

p = p- Pm: r = r/rm , & = log = - log . ,1 ° . (3.30) 



,'a' 4 4q;'2 

Next we proceed by solving numerically the equation of motion for the classical embedding of the 
probe x{p) expanding at p = = log(r*) 



sin = - + - + O ( ^ ) with m = mo/rm & c = c/r^; , (3.31) 



where as we discussed above mo & c are related to the bare mass and the fundamental condensate 
of the Af — 2 hypermultiplet corresponding to the probe brane. Next we vary the parameter e* and 
generate plots of — c versus rh, which we present in figure 1. 

The first plot corresponds to a very small value of e*, namely 10~^. One may expect that for such 
a small value the qualitative and quantitative behavior of the system would be the same as for 
the case of pure AdSs x space-time, studied in [3]. Indeed our results confirm that. The other 
plots in figure 2 correspond to = 0.25, 0.5, 0.75. One can observe that qualitatively the theory 
has the same properties. For sufficiently low bare mass (strong magnetic field) there are multiple 
phases forming a spiral structure near the origin of the — c versus m plane. Only the lowest positive 
(m > 0) branch of the spiral corresponds to a stable phase [22]. Furthermore at vanishing bare 
mass the theory has a non- vanishing negative condensate which spontaneously breaks a global U{1) 
R-charge symmetry. As one can see the only effect of the backreaction is to lower the value of this 
condensate. 

In figure 2 we have presented a plot of the symmetry breaking condensate (at vanishing bare mass) 
as a function of the parameter e*. We have normalized the plots in units of the condensate at 
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vanishing e*. Note that one can vary e* either at fixed (fixed magnetic field) or at fixed infrared 
separation Lq = rmin (the radial distance above the origin at which the wrapped by the probe 
brane shrinks). As one can see in both cases the fundamental condensate decreases as we increase 
e*. Note that in figure 2 we keep the value between zero and one, even if it is possible to extend 
this range. The reason behind this choice is the separation of scales, which keeps the finite cut off 
well below the Landau pole only in the limit -C 1. 

4. Probing the fully backreacted background 

In this section we introduce a probe D7-brane to the perturbative background obtained in section 2. 
On the field theory side this corresponds to an introduction of additional fundamental fiavors. The 
fundamental fields couple to the background i?-field, which on the field theory side corresponds to 
an external magnetic field. Therefore, in the limit of vanishing bare mass we expect to observe the 
phenomenon of magnetic catalysis of mass generation according to [3]. The novel feature of our set 
up is the presence of smeared backreacted fiavors which number is controlled by the parameter e*. 
Our immediate goal is to study the effect of the presence of the sea of massless fundamental fields 
on the process of mass generation. Note that in this set up the external magnetic field couples to 
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fixed r„ 



fixed Lo 




Figure 2: Plots of the dependence of the scaled fundamental condensate {qq) e, / {qq) o versus the small 
parameter e^,. The left plot corresponds to varying e^: at fixed (fixed magnetic field) and the right plot 
corresponds to fixed infrared separation Lq (fixed constituent mass). 



the smeared massless fundamental flavors. This is to be contrasted to the case studied in section 
3.2, when the external magnetic field couples only to the flavor fields introduced by the probe brane. 
Our investigation indicates a significant change in the physical properties of the system. 

4.1 Classical embeddings 

Let us begin by considering a probe D7-brane extended along the t,x\a,9,ip,(f> coordinates and 
having a non trivial profile along the coordinate x- The effective action of the probe in the Einstein 
frame is 

5 = -Tt y rfV^-det(^ + e-*/252)+ Tr J P[C^s) - 5(2) A Qg)] , (4.1) 

where C(8) and C(6) are background Ramond-Ramond forms sourced by the smeared flavor branes. 
One can show that the effective lagrangian is 

£(°) oc le^'S'F^ cos^ ¥ cos' ? + 4 sin^ ^) 1+ ^ A ■ ^ ^ ^ 



8 2V 2 F2 2j \ Ah^S^F'^ J \ b"^ 

+ |e-6^5«(l + ^^)cos^|. (4.2) 

We find it convenient to introduce the radial variable r related to a via (2.25). In order to obtain 
numerical solutions for the classical embedding, we define as usual the following dimensionless 
variable 

r = — , with Tm = T . (4.3) 

It is convenient to write the dimensionless action in terms of the following auxiliary functions 

= -/i(r) (cos^ I + f2{rf sin^ |) ' ^/l + Mr^x" + W) cos' | , (4.4) 
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where 



1 ^ 







)■ 


df 



h{f) = (2bS'F 



da 



df 



-1 



' ^ ' 32 I f%'^ J df 



(4.5) 



and ^ ,b , S , F , H ,a are related to the functions defined in equations (2.25)-(2.32) via 

<|)(f) = $(rT„) - $^ , 6(f) = 6(fr^) , S{f) = S{fr^)/r^ , 
F = F{fr^)/rm , ^(f) = H{fr^)/H^ , a(f) = cr(fr^)r^ . 



(4.6) 



One can show that a smooth sohition to the equations of motion derived from (4.4) vahd near the 
point that the S'^-sphere, wrapped by the probe brane, shrinks {xiTmin) — tt) is 



x(0 



TT 



r — Tr, 







a — 



8/1/2 



h{-2h + /i/s/^ + f2{hf[ + /i/^)) 



(4.7) 



The approximate analytic solution (4.7) can be used to fix the boundary conditions for a numerical 
shooting technique at f = fmin + ^ for small values of 5. 

A crucial point is that due to the infrared divergency of the perturbative solution, (2.25)-(2.32), 
this is not valid below a certain radial distance fiR^e^) above the origin (f = 0). The probe reflects 

that by having its tension vanishing at f7R(e*) but in fact it is the Jacobian |? which is vanishing 
at this value of the radial coordinate. As one may expect if we increase the number of backreacted 
flavors the radius f//j(e) grows. Therefore we start loosing some of the probe brane embeddings 
corresponding to low bare masses. Fortunately this effect is somewhat softened by the effect of 
magnetic catalysis of mass generation which leads to infrared separation of the probe branes (the 
probe closes at some distance above the horizon). To visualize this behavior in flgure 3 we have 
presented plots of D7-brane embeddings in (r, x) polar coordinates for a range of different infrared 
separations (Lq = fmin) and different number of backreacted flavors (e* oc A* -jf). 
The flrst plot corresponds to the absence of backreacted matter studied in [3]. As one can see there 
are two classes of embeddings: those which asymptote to negative (positive) separation at large f. 
One can show that [3] only the second class corresponds to stable embeddings. The two classes are 
separated by a critical embedding (depicted with a red curve in flgure 3), which has zero separation 
at large f (vanishing bare mass) and flnitc separation in the infrared (dynamically generated mass). 
The second and the third plots correspond to small (non- vanishing) values of the parameter e*. As 
one can see due to the infrared divergencies some of the embeddings with small separation in the 
infrared are lost. However the critical embedding is still present and we can study how its separation 
in the infrared (related to the dynamically generated mass) depends on the number of backreacted 
flavors. 
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Figure 3: Plots for the profiles of D7-brane embeddings having different infrared separations for vari- 
ous values of the perturbative parameter e*. The red curves correspond to embeddings with vanishing 
separation at r = r*, which we interpret as describing fundamental fields with vanishing bare mass and 
non-vanishing constituent mass. 



The last plot in figure 3 corresponds to a sufficiently large number of e*, at which the radius of the 
region which is not accessible by the probe brane has increased so much that the embedding with 
the lowest separation in the infrared asymptotes to positive separation at large f. 
At first sight it seems that the validity of the perturbative background is not sufficient to detect 
a qualitative change of the theory as a function of the parameter e^. Furthermore unlike the case 
of the supersymmetric background studied in section 3, we do not have an analogous holographic 
renormalization procedure. This is why we cannot justify the validity of the AdS/CFT dictionary 
and in particular (3.27). 

Fortunately one can use the properties of the D7-brane embeddings in the infrared, such as the 
infrared separation Lq to study the constituent mass of the fundamental fields Mq. Indeed a detailed 
study of the behavior of the infrared separation I/o(e*) of the critical embedding as a function of 
the parameter e* is presented in figure 4. 

The blue curve corresponds to the function Lo^e^). The dashed growing curve corresponds to 
the radial distance fm^e^) below which the perturbative solution for the background geometry 
cannot be trusted. The study indicates that for sufficiently large values of e^. the function I/o(e*) 
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Figure 4: A plot of the infrared separation Lq (at zero bare mass) versus oc A* ^ 



becomes multivalued through a diverging slope at Lq k, 0.326. This in turn implies that the critical 
embedding becomes unstable at this point. The study of the meson spectrum supports that. 
It is plausible to interpret the instability of the probe D7-brane as reflecting an instability of 
the background. Indeed if we consider a single "fiducial" embedding, representative for the smeared 
massless flavor branes such an embedding would be unstable in the probe limit due to its coupling 
to the background 5- field [22]. Therefore we would expect that the background constructed by 
smearing many such embeddings would be unstable. 

Then a natural question arises. Is there a critical value of the B-field at which the instability is 
triggered or the background is unstable at any non-vanishing magnetic field? In the probe limit at 
zero bare mass the only physical scale is the scale corresponding to the non- vanishing 5-field (i?*) 
therefore one would expect that the background is unstable at any H^, > 0. On the other side in the 
backreacted case the dilaton field is running and the theory develops a Landau pole. This suggests 
that even at vanishing magnetic field the backreacted theory has a physical scale associated to the 
finite UV cut off needed to keep the relevant energy scales below the Landau pole. Therefore one 
would expect that there is a finite Her > at which the background becomes unstable. 
Let us consider again the physical process described in figure 4. Qualitatively the quantity Lq can 



^°Note that unlike the critical embedding, the smeared massless flavor branes correspond to trivial (x = 0) 
embeddings, which are unstable in the presence of non-vanishing _B-field 
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be thought of as proportional to the ratio of the dynamically generated mass of the fundamental 
fields introduced by the probe D7-brane Mq and the energy scale corresponding to the external 
magnetic field ^/H^, namely Lq oc Mg/ y^TT^. At any fixed value of we have three independent 
energy scales: Mg, H^, & Her- In order to compare states at different values of e^, we will keep the 
dynamically generated mass of the quarks fixed. The data in figure 4 suggests that as we increase 
the value of Lq decreases. At fixed value of Mq this is equivalent of increasing the external magnetic 
field if*. Assuming that the critical value Her varies slowly enough with e* we conclude that at 
sufficiently large values of > e^r the external magnetic field H^, is H^, > H„. and the gravitational 
background is unstable. It is plausible to interpret the instability of the probe D7-brane associated 
to the diverging slope of I/o(e*) in figure 4 as refiecting an instability of the gravitational background 
for > Her- 

Note that the limited validity of our perturbative solution in the deep IR prevent us from directly 
studying the stability of the background. For example by studying the spectrum of minimally 
coupled scalar field. It would be interesting to address this problem at finite temperature, where 
the IR behavior of the supergravity solution is regularized by the non-extremal horizon. We leave 
such studies for future work. 

The instability of the probe D7-brane could be the sign of a "new" physical effect, called "super- 
conducting vacuum"^^. According to [23], the quantum vacuum (i.e., an empty space) may become 
a superconductor (in the usual electromagnetic sense) under a strong enough magnetic field. This 
magnetic field forces the electrically charged vector mesons to condense via a tachyonic instabil- 
ity, which in turn implies electromagnetic superconductivity. Calculations in a variety of models, 
supporting this idea, can be found in [24, 25, 26]. 

4.2 Meson spectrum 

In this section we analyze part of the spectrum of meson like excitations of the dual gauge theory. 
To this end we perform a semi-classical quantization of the probe D7-brane embcddings studied in 
section 4.1. In order to obtain the spectrum of the corresponding quantum fiuctuations we expand 
systematically the effective action of the probe D7-brane to second order in a'. 

4.2.1 Quadratic fiuctuations. 

For a detailed study of the light meson spectrum of the theory we need to consider the quadratic 
fiuctuations of a D7-brane and study the corresponding normal modes, [11]. The relevant pieces of 
the action are 



7-- J rf'x e*^- det [Gab + e-? {Bab + Fab)] + J P[C(8)-F(2)AC(6)+^F(2) A F(2) A C(4)] , (4.8) 
where C(4) , C(6) & C(s) are defined by: 

C(4) = C'(4) - C(2) A S(2) , C(6) = C(6) - -6(2) A C(4) , C(8) = C(8) - -B(2) A C(6) (4.9) 

However we should point out that though plausible our studies based only on a probe brane calculation alone 
are not conclusive. 

^^We thank Maxim Chernodub for pointing out to us this possible interpretation. 
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and C(4) is the magnetic dual of the RR potential C(4) 

(7(4) = 

We will consider fluctuations of the form 



1 .V 

— Qc sin 9 cos ^ dO A d(p A dip A dr , 



X = Xo(<7) + 2W5x(r) & T = 2W5r(r), 



(4.10) 



(4.11) 



where the indices a, & = 0, 1, . . . , 7 run along the worldvolume of the D7-brane and expand (4.8) 
to second order in ol . Before presenting the relevant terms in this expansion it is convenient to 
introduce S &i J, which are symmetric & antisymmetric matrices respectively, in the following way 



\Kb\\-' = s + J. 



(4.12) 



The non-zero elements of those matrices are 

_ ^11 _ Q-l ^22 _ ^33 _ 



G 



22 



C<Ta 



G 



f~< r< r~<2 



r' c' 



ri-i ode _ ^-1 



r' r' c"^ ' 



e-2H 



Gl^ + e-*ii'2 



(4.13) 
(4.14) 



with 



Gn = gf^ , G22 = 9^^ , G,, = + xM' , 

/-< _ (0) ri 
^ee — 900 , ^ipiii - 

The second order terms in a' expansion of the action (4.8) are 



(0) r - n(°) r - n(0) 



(4.15) 



T7 



(0) _ J0)2q^^ _Ort(0)„(0) _ {0)2 



9iPt 



S'^'dadTd.dr , 



S'^'dJxd.Sx + - fia) 5x' , 



-4xSr = TrJ'^SxdaSr , = ^ v^'5""''5'''^a(.i^a';,' + Y ^[^4] e"^""^ , (4.16) 



32 



cos 



StFu . 
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where the indices m, n, o, p run in the transverse directions of the D3-branes. The auxihary functions 
f,g^ J"' 1 appearing in the above equations, are defined as follows 



(4.17) 



(0) 



where C^^\d^j^ are the complete classical lagrangian of the probe brane (4.2) and the Cq contri- 
bution of the Wess-Zumino part respectively. Note that the function J"' is non-vanishing only for 

a = , 

Note that if one restricts the fluctuating modes to depend only on the time and the holographic 
directions the coupling term in the effective lagrangian vanish and we can consistently focus only 
on the fluctuations along x- 



4.2.2 Fluctuations along x 

We consider the ansatz 



(4.18) 



The only relevant part of the effective action is the term C,f^^-^- We find it convenient to use the 
radial coordinate ria) related to a via equation (2.25). Using the dimensionless notation defined in 
equations (4.3) and (4.6) we obtain the equation of motion for 77(f) 



d, 



''-DBliS 
(I + /3VT 



''-DBli5 ~2 



1 + fix" 



UJ — 



/92£(o) _ o / fix' a m 



77 = 0, 



where 



UJ = 



4^2 



2r 



(4.19) 



(4.20) 



In order to solve numerically (4.19) we need to impose proper boundary conditions at f = rami = Lq. 
Using (4.7) for the classical embedding one obtains the following asymptotic form of (4.19) 



■q"{r) + - — — 
The general solution of equation (4.21) is 

r]{f) = — 



-r]'{f) + 



4 {t T"min) 



Mr) = 



(4.21) 



Co 



(4.22) 
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Requiring renormalizability for 77(f) 



J dfV^W) ri{f ) \ 00 , (4.23) 



and the fact that v ^ ^'■"^ ~ (^^ — ^mm) on can see that C2 should vanish. This suggests the following 
substitution in equation (4.19) 

V(r) = -^ML= , (4.24) 



r — 



where the function ({r) is analytic near Trnin- The resulting equation of motion for ({f) is of the 
form 

\ I ' mm / \ ' ' mm ^ V ' min ) J 

where Cq ,Ci ^ B are functions of f and x(f) regular at f = f„ii„. Equation (4.25) has predetermined 
boundary conditions at Tmin- Therefore one can expand C(f) near f„jj„ and solve equation (4.25) 
order by order. The obtained approximate solution for ({r) can be used to fix the boundary 
conditions for a numerical shooting technique. The resulting numerical solution depends on the 
parameter a). Imposing Dirichlet boundary condition at f = f* quantizes the spectrum of a). 
We used the procedure outlined above to generate plots of the first three excited states of the 
spectrum as a function of the parameter e^,. The resulting plots are presented in figure 5. One 
can see that the ground state becomes tachyonic exactly at the critical value of ^ 0.058 where 
the slope of Lq versus in figure 4 diverges. This confirms that at this point the probe branes 
become unstable. As we commented in subsection 4.1 we interpret this instability as reflecting an 
instability of the background. This suggests that the dual gauge theory is unstable for sufficiently 
strong magnetic fields. This is to be expected because the backreacted massless fiavors are in a 
phase with vanishing constituent mass which is disfavored by the external magnetic field. 



5. Conclusions 

In this paper we propose a string theory dual to a 1+3 SU{Nc) N — ^ SYM theory coupled to 
Nf massless fundamental flavors in an external magnetic fleld. Our motivation is to undertake the 
first steps towards an unquenched holographic description of magnetic catalysis of mass generation. 
By construction, our background corresponds to a phase of the fiavored theory with vanishing con- 
stituent mass. This phase is disfavored by the external magnetic field, therefore one would expect 
that for sufficiently strong S-field the supergravity background would become unstable. Unfortu- 
nately the infrared singularity of our perturbative solution prevents us from directly studying the 
stability of the background by studying the spectrum of its quasi-normal modes. To circumnavigate 
this limitation we study the properties of an additional probe D7-brane. 
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Figure 5: Plots of the first three excited states of oj versus oc Nf/Nc- One can see that at the critical 
value e* w 0.058 the ground state is tachyonic. 

In Section 3.1 we study the properties of a supersymmetric probe D7-brane in the limit of vanishing 
-B-field. In this way our background reduces to the supersymmetric one obtained in [15]. We 
consider a holographic renormalization of the probe brane "on-shell" action in the spirit of [21], but 
in the case with backreacted flavors. Our studies reveal a remarkable factorization of the dependence 
of the "on-shell" action on the perturbative parameter counting the number of backreacted flavors 
(e*). This factorization suggests that the holographic renormalization performed in [21] can be 
implemented, at least formally, for the background [15]. Next, a systematic expansion in e*, for 
e^, ^ 1 & keeping the finite cut off of the theory sufficiently far bellow the Landau pole, provides 
a regime of validity of the renormalization procedure. This suggests that the usual AdS/CFT 
dictionary holds at the UV scale fixed by the finite cut off. 

In Section 3.2 we move one step forward our study and introduce a non-supersymmetric probe D7- 
brane to the supersymmetric background of [15]. This probe brane has a fixed U{1) worldvolume 
gauge field corresponding to an external magnetic field coupled to the fundamental fields introduced 
by the probe brane only. Using the AdS/CFT dictionary proposed in the previous section we study 
the effect of mass generation and its dependence on the number of backreacted fundamental flavors. 
Qualitatively the physical picture remains the same compared to the one without backreacted 
flavors, [3]. We contrast these results to those of a magnetic flled coupling to all fundamental 



26 



degrees of freedom. 

Finally in Section 4 we consider a D7-brane probe in the supergravity background obtained in 
Section 2. Now the i?-field, corresponding to an external magnetic field, couples to both the probe 
and the backreacted fundamental degrees of freedom. The study of the meson spectrum shows that 
for sufficiently strong magnetic field the probe brane becomes unstable. It is plausible to interpret 
this instability as reflecting an instability of the supergravity background. As commented above 
this is anticipated since the external magnetic field disfavors a phase with vanishing constituent 
mass of the fundamental fields. We speculate that the stable phase at vanishing bare mass would 
correspond to a supergravity background obtained by smearing "fiducial" embeddings with a non- 
trivial profile along the radial coordinate having finite separation in the infrared corresponding to 
dynamically generated constituent mass. The construction of such a background is one of the main 
directions for future studies that we intend to pursue. 

A possible way to improve the IR properties of the gravity dual constructed in Section 2 is to 
consider the non-extremal case, when the zeroth order solution is the AdSs x black hole. In this 
case it would be possible to study directly the spectrum of quasi-normal modes of the background 
and verify the anticipated instability for sufficiently strong magnetic field. One could also introduce 
an additional probe brane and perform a study analogous to the one considered in Section 4.2 of 
the present work. Naturally, one would expect that such a probe would become unstable when the 
energy scale of the magnetic field is sufficiently larger than the energy scale of the finite temperature. 
One can then compare the critical parameters obtained from studying the spectrum of quasi-normal 
modes of the background and from studying fluctuations of the probe brane. Such a study could 
provide an indirect check of the interpretation of the instability of the probe brane cosidered in 
Section 4.2.^*^ We leave such studies for a future work. 
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A. Equations of motion 

The equations of motion produced from (2.1) are [27] 



2 3! 



2 

(3) 



1 1 

+ 23!^" 



2 

(3) 



8 

1 1 (5(5' 

d[^d$] = e=^*F(i) A*F(i) + -e*F(3) A*F(3) - -e-" H^s^ A^H^s) - 2^ 



(3) - ^ ^2 , 



d [e* * F(3)] = -//(3) A F(5) + - A ^2 , 

d [*F(5)] = dF(5) = i^(3) A F(3) - ^ ^2 , 

d [e-* * i/(3)] = e*F(i) A ★F(3) - F(5) A F(3) 



ST 



(A.l) 

(A.2) 
(A.3) 
(A.4) 
(A.5) 

(A.6) 



where the last term in (A.6) is an eight form and denotes the derivative of the smeared DBI action 
with respect to T 



d^a;e*A/-det(G + e-*/2^) ^ / d^"a;e*./- det(G' + e-*/2j^)|02 



(A.7) 



The symbol IQ2I that appears in (A.7) denotes the modulus of the smearing form and has the 
following expression 



in-, 



2Qf 



The Bianchi identities for all the forms of the background are 

dF(3) = //(3) A F(i) - ^5 J" A Q2 , 
d//(3) = 0. 

The Bianchi identities and the equations of motion for the worldvolume gauge fields are 

6Sfi 



0. 



(A.8) 

(A.9) 
(A.IO) 
(A.11) 

(A.12) 
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Plugging the ansatz of (2.7) in the Bianchi identities (A. 9), (A. 10) and (A. 11) as well as the 
equations of motion (A. 3) and (A. 5), trivially satisfies them all. The equations of motion for F3 & 
H^, namely (A.4) & (A.6), will give us (2.13) & (2.21) respectively. 

The equation of motion for the dilaton, namely (A. 2), will give us (2.20). The contribution from 
the last term in (A. 2) is the following 



6S 



fi 



Ae^Qf 1 + 



262 



1 + 



(A.13) 
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The Einstein equations, namely (A.l), will give us (2.16), (2.17), (2.18) & (2.19) together with the 
constraint (2.22). This last one is coming from the 44 (frame) component of the Einstein equations. 
The contribution from Tfi that appears in (A.l) is the following 



2^2 SS 



-9 5g- 



9s e 



_0(2)o(2)p 



(A.14) 



with 



= V,. + e-? ^ {5i^5l - 5^^5^,) , (/., = 0, 1, . . . , 9) 



and S is the symmetric part of the inverse of £ 

1 



diag{-l,l. 



JJ2 fl J 1) 1) 1) 1) 1) ^} ■ 



(A.15) 



(A.16) 



Finally for the worldvolume gauge fields, we have the first of the equations in (A. 12) trivially satisfied 
while the second splits into 



d 



+ d 



6S^ 







(A.17) 



The first term in (A.17) can be easily evaluated if one differentiates (A.6) and plugs in the ansatz 
(2.7). One can show that it is of order d [0{dA)]. The second term is proportional to the expression 



d [C(6) A + 0{dA)] d [0{dA)] 



(A. 18) 



The bottom line is that the whole expression in (A. 12) is of order d[0{dA)] and hence one can 
consistently set the D-brane worldvolume gauge field to zero (^4 = 0). 



B. Technical details 



In this appendix we will present some technical details on the computations of Section 3. 
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B.l Exact On-Shell Action 

The "on-shell" action —Sci/{Ma''^) of a supersymmetric embedding witli bare mass parameter e^'^ 
in the supergravity background described by equations (3.17)-(3.20) is given by 



/ 



Pq 



dp 



e<f*(e2p - e''P'')[Qe,{e''P - e''P''){l + e*(| + - p)) + 4(1 + e,(p, - p)){<oe^P{l + e,{p, - p) + e.e^P" 



32 X 6(1 + + p,- p))^(l + e.(p* - p))^ 



(B.l) 

The definite integral in (B.l) can be solved in a closed form using the relatively simple dependence 
on the parameter e''«. To this end one defines rriq = e^'' and takes the second derivative with respect 
to rrig = e^f" 



d'^Sci[p*,m, 
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Furthermore equation (B.l) suggests that Sci\p =0 and one can verify that dirn^ySd 
Therefore we can integrate (B.2) and obtain the following expression for the classical action 



P*=Pq 



Afa'^ 32 



(l + |)^(e^--e^^^)^ 



(B.2) 
= 0. 

(B.3) 



which is exactly (3.23). 
B.2 The Function h 



In order to evaluate the effective action from (3.29) we need to integrate the equation of motion for 
h{p) in (3.20). Changing variables to dimensionless ones along (3.30), the equation of motion for 
h{p) becomes 

dh 



dp 



1 + e* ( - + p* - p 



(B.4) 



A short integral expression can be given in terms ofa; = l/e + l/6 + p* — p 



X3 



2 

e*3 



4.x \ I 

2 I with X* = x{p^) = h - 



6 



(B.5) 



The constant of integration in (B.5) is fixed by requiring h{p^) = e ^''*. Using that 

„4x 



dx- 



■2n 

e 3 



r|i-4. -r 1 



+ const , 



X3 23 

where r(a) and r(a, z) are the complete and incomplete gamma functions^^, we have 



1 + 



4g-|-|{i-i.) 



(2e. 



^,142,, A ^ /I 4 2 
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Note that the right-hand side of equation (B.6) is real for a; > 0. 



(B.6) 



(B.7) 



30 



Note that for p G (— oo,p*) & > h{p) is real. 
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